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We have carried out nuclear quadrupole resonance (NQR) measurements in the filled 
skutterudite compounds Pr(Osi_ a; Ru a ;)4Sbi2 (x = 0.1, 0.2), in order to gain insights into the 
symmetry of the superconducting gap function. Upon replacing Os with Ru, the spin-lattice 
relaxation rate 1/Xi becomes proportional to temperature (T) at low T far below T c , and 
the magnitude of (\/TiT)i ow ^t increases with increasing Ru content. These results indicate 
that a finite density of states is induced at the Fermi level by the impurity, and thus suggest 
that there exist nodes in the gap function of PrOs4Sbi2 . 
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The superconductivity at T c =1.85 K discovered in the filled skutterudite compound 
PrOs4Sbi2 has attracted much attention. 1 ' 2 This is the first Pr-based heavy fermion su- 
perconductor, in which the heavy mass has been suggested by the large specific heat jump 
AC/T C ~500 mJ/K 2 - mol at T c , 1,2 and directly confirmed by de Haas-van Alphen (dHvA) 
experiments. 3 The crystal electric field (CEF) ground state is a Ti singlet, which is separated 
by the first excited state of the 1^ triplet by a gap of Acef=W K. 4 ~ 9 Because of this small 
A-CEFi the role of quadrupole moment fluctuations arising from the T 4 state has attracted 
much attention. It has been speculated that the interaction between the quadrupolar moments 
and the conduction-electron charges may be responsible for the heavy mass, thus representing 
a charge-scattering version of the Kondo effect. 10 The contrasting behavior that the isostruc- 
tural compound PrRu4Sbi2 with Acef=70 K has a much smaller m* 11 ' 12 has increased such 
expectation. 

One of the first steps toward understanding the Pr-based heavy fermions may be the 
determination of the gap symmetry of the superconductivity. Previous nuclear quadrupole 
resonance (NQR) measurements have revealed the unconventional nature of the supercon- 
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ductivity. 13 No coherence peak was found in the spin-lattice relaxation rate 1/T\ just below 
T c , which suggests non-BCS-type superconductivity and is in contrast to the s-wave bahavior 
subsequently observed in PrRu4Sbi2- 14 At low temperatures, 1/T\ does not follow a power- 
law T-dependence, which is in contrast to the observation of a T 3 variation in other known 
heavy fermion superconductors. 15 Together with the exponential decrease in the penetra- 
tion depth found by muon spin relaxation, 16 it has been suggested that the superconducting 
gap is isotropic. However, oscillations with respect to the magnetic field angle have been 
found in angle-resolved thermal conductivity measurements, which suggests that the gap is 
anisotropic. 17 

In this study, we use an impurity as a probe for the gap function. We replace Os with Ru, 
and study its effect on the superconducting density of states (DOS) using the NQR technique. 
It has been known for some time that a non-magnetic impurity is a powerful "smoking gun" 
for the gap function. In s-wave superconductors, a paramagnetic impurity produces a bound 
state within the gap (Yu-Shiba state), 18 ' 19 and non-magnetic impurities smear out the gap 
anisotropy, if any. 20 In contrast, if there are nodes in the gap function, a finite DOS is brought 
about at the Fermi level, 21,22 which can be detected experimentally. 15 

Single crystals of Pr(Osi_a;Ru x )4Sbi2, (x=0.1 and 0.2) were grown by the Sb-flux method. 
T c 's are 1.6 K and 1.4 K for x=0.1 and 0.2, respectively, in agreement with those previously 
reported. 23,24 For NQR measurements, the single crystals were powdered in order to allow the 
radio-frequency magnetic field to penetrate into the sample. The NQR measurements were 
performed by using a phase coherent spectrometer. The spin-lattice relaxation rate, 1/T±, for 
123 Sb were measured using the saturation-recovery method. Data below 1.4 K were collected 
using a 3 He/ 4 He dilution refrigerator. 

Figure 1 shows the NQR spectra for the x=0.1 and 0.2 samples. For comparison, the 
data for the two end-member compounds 13, 14 are also shown in Figs 1(a) and 1(d). In the 
Ru-doped samples, in addition to the spectrum peaks that are located at the same positions 
as those for pure PrOs4Sbi2, there appear new sets of peaks that are characterized by uq=26 
MHz and the asymmetry parameter 77=0.43 for 123 Sb. Here, vq and r/ are defined as uq = 
Vz = 2J(2J-i)fe e2( ^^^ an< ^ ^ = \ Vx ~ v y\l v z with Q being the nuclear quadrupolar moment 
and 0(a = x, y, z) being the electric field gradient at the position of the nucleus. 25 It should 
be emphasized that these peaks (labelled as "peaks 2" in Table I) are different from those 
for pure PrRu4Sbi2- These new peaks are due to the alloying effect, which indicates that the 
alloying is homogenous and there is no phase separation. Such a feature has been seen before 
in other heavy fermion compounds. For example, in Ir-doped CeRhlns, upon replacing Rh 
with Ir, there emerges a new set of NQR peaks that is different from that for CeRhlns or 
Celrlns. 26 Nonetheless, 1/T\ measured at the new set of peaks shows the same T-dependence 
as that measured at those corresponding to CeRhlns or Celrlns, indicating that the new set 
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Fig. 1. Sb-NQR spectra of Pr(Osi-a;Rua;)4Sbi2 at 4.2 K. The data for x = and x = 1 are from 
Kotcgawa et al. 13 and Yogi et al. 14 In (b) and (c), the spectra indicated by the thin solid lines 
(broken lines) correspond to 123 Sb ( 121 Sb) signals seen in PrOs4Sbi2, while the spectra indicated 
by the thick markers are the new sets arising from Ru doping. 

of peaks is not due to phase-separated regions but rather to the homogenous alloying. 

The nuclear spin-lattice relaxation rate 1/T\ was measured at the peaks at 49.9 MHz 

(±3/2 < ► ±5/2 transition) and 42.2 MHz (±1/2 < ► ±3/2 transition). Both measurements 

yield the same results. Even for x=0.2, no influence on T\ due to the overlapping of the 121 Sb 
line was seen in the measurements at the 49.9 MHz peak, probably because the Hi we used for 
observing the echo is much smaller than the frequency difference between the 121 Sb and 123 Sb 
peaks, which is about 0.8 MHz, and the 121 Sb peak has a small full width at half maximum 
(FWHM) of 0.5 MHz. Above T ~ 100 K, the nuclear magnetization can be fitted excellently 
using the expected theoretical rate equation. 27 However, below T ~ 100 K, it cannot be fitted 
using the theoretical curve with a single component of T\. This is also true at the peak at 42.2 
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Table I. NQR frequency vq and asymmetric parameter 77 for 123 Sb in Pr(Osi_ x Ru x )4Sbi2- Also 
shown are the CEF gap A CEF and the heavy quasiparticle enhancement factor (see text). 
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Fig. 2. Nuclear magnetization recovery curve. The broken and dot-dashed curves indicate the recov- 
ery curves due to the fast and slow relaxation components, respectively. The solid curve indicates 
the sum of the two components. 



MHz that is free from overlapping by other transition lines. Figure 2 shows the decay curve 
of the nuclear magnetization at T=6 K. As will become clearer later, the inhomogeneous T\ 
is an intrinsic property of the alloyed samples. The one-component behavior for T > 100 K is 
rather a coincidence, since T\ is the same for both PrOs4Sbi2 and PrRu4Sbi2- 

We then attempt to fit the low-T nuclear magnetization with two T\ components, namely, 

M - M(t) _ 
Mo ~ 

£i=S,L°i [0.0762exp (-f|) + 0.0146 exp (-^f) + 0.9092 exp (-^f )] ,(1) 

where Tf means the short component due to fast relaxation, and Tf, the long component 
due to slow relaxation. The solid curve in Fig. 2 shows such fitting for x=0.1, with 05=0.8 
and aL=0.2 that do not depend on temperature appreciably. 

We first discuss the main component, l/Tf, whose temperature dependence is shown in 
Fig. 3. Compared to the results for the end-member compounds, several trends can be seen. 
First, above T ~ 100 K, there is only one component of T\ (namely, Tf = Xf), and the data 
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Fig. 3. Temperature dependence of l/Tf in PrtOsi-^Ru^Sb^ (a; = 0.1, 0.2). Data for PrOs 4 Sbi 2 13 
and PrRu 4 Sb 12 , which is a BCS superconductor, 14 are also shown for comparison. The straight 
lines indicate the TiT=const relations. 



for all x merge into a same line. This suggests that the relaxation at high T is not governed by 
the electronic state that is sensitive to the transition metal element, 28 but due to other degree 
of freedom. Second, below T ~ 100 K, l/Tf is close to that for pure PrOs4Sbi2, but decreases 
as x increases. Since the slow component l/Tf is very close to that for pure PrRu4Sbi2, as will 
be shown later, we consider that l/Tf origenates from the Sb sites that are located far from 
Ru. Third, and most importantly, below T ~ 1 K, l/Tf is proportional to T; the magnitude 
of 1/TfT increases with increasing x. 

The l/T 1 s 'r=const. relation indicates that a finite DOS is induced by the impurity. 1/T\ 
in the superconducting state may be expressed as 



T 1 (T = T C ) 



+ M s (E) 2 )f(E)(l-f(E))dE, 



(2) 



Ti k B T c 

where N S (E) = NqE/(E 2 - A 2 ) 1 / 2 is the superconducting DOS with A being the super- 
conducting gap, M S (E) = NqA/(E 2 — A 2 ) 1 / 2 is the anomalous DOS due to the coherence 
factor, 29 iVo is the DOS in the normal state and f{E) is the Fermi function. In the T=0 limit, 
1/Ti is dominated by iV s (0) 2 . The present results indicate that N S {Q) is finite in the Ru-doped 
samples and increases with increasing x. This result cannot be accounted for by an s-wave 
gap. 18-20 It is emphasized that only if there is a sign change in the gap function, namely, there 
are nodes in the gap function, can the impurity induce the finite DOS. 20 Thus, our results 
provide strong evidence for the existence of nodes in the gap function. The seemingly expo- 
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Fig. 4. Suppression of T c as a function of the residual DOS induced by the impurity, where T c q is the 
transition temperature in the absence of impurity scattering. The curves indicate the calculated 
results obtained by Miyake. 21 ' 22 




nential decrease in \/T\ with respect to T in pure PrOs4St>i2 may be an effect obscured by 
the double superconducting transitions. 30,31 Let us now estimate the impurity- induced iV s (0), 
denoted as N res hereafter. The values of N res calculated from the relation of 

Nres = I (T(T) Tc 

N V (T(T) low _ T 1 ' 

are shown in Fig. 4. The curves indicate the calculated results obtained by Miyake for the gaps 
with line-nodes (axial) and point-nodes (polar). 21 ' 22 The experimental results agree qualita- 
tively with the theoretical results. That the experimental data fall below the theoretical curve 
is probably due to the depression of the pairing force by the substitution of Ru, which is not 
included in the theoretical results. We will discuss this point later. Although we are unable to 
distinguish between line nodes and point nodes in Fig. 4, point nodes nonetheless seem more 
plausible, in view of the relatively weaker suppression of T c by the impurity with respect to 
x. 23 

For completeness, we show the temperature dependence of in Fig. 5. Above T c , l/Tf 
shows a temperature variation very close to that of PrRu4Sbi2- This 1/T± can be assigned to 
come from the Sb sites that are located close to Ru. Below T c , 1/T[ decreases rapidly, with 
no coherence peak. 

Finally, we discuss the high-temperature behavior oil/Tf . Figure 6 shows the temperature 
variation of the quantity 1/T\T. Below T ~ 6 K, 1/T\T decreases, leaving a peak at around 
T ~ 6 K. Such a decrease becomes mild as x increases. Since the CEF gap Aqef is small, 
Kotegawa et al. have analyzed the low-T data by decomposing 1/T]T into two parts, 13 namely, 

(2) 

the contribution due to the excitation to the T 4 state, and that due to heavy quasiparticles: 
\jT\T = axexp(-A C EF/k B T)+l3x0.7 s^K- 1 , where 0.7 s^K- 1 is the 1/TiTfor LaOs 4 Sbi 2 
and the parameter y//3 represents the mass enhancement factor of the heavy quasiparticles. 
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Fig. 5. Temperature dependence of the minor, slow relaxation component, 1/Tf. The straight line 
indicates the l/7iT=0.7 s K relation found in LaOs4Sbi2- 



The application of the same analysis of Kotegawa et al. to the x=0.1 and 0.2 samples, as 
shown in the inset of Fig. 6, shows that Aqef increases with increasing x, while f3 decreases 
with increasing x. The results are summarized in Table I, which are in good agreement with 
those reported by Frederick et al. who deduced Aqef from the susceptibility and the mass 
enhancement factor from the specific heat coefficient 7. 23 As has been discussed by several 
authors, 14,23 the increase in Aqef is responsible, at least partly, for the decrease in T c in 
going from PrOs4Sbi2 to PrRu4Sbi2; thus, the discrepancy between the theoretical curve and 
the experimental data seen in Fig. 4 can be ascribed to the increase in Aqef- 

In conclusion, we have studied the impurity effect on the superconductivity in the filled 
skutterudite heavy fermion superconductor PrOs4Sbi2 using the NQR technique. We find that 
replacing Os with Ru brings about a finite density of states at the Fermi level, which increases 
with increasing Ru content. Our results provide strong evidence for the existence of nodes in 
the gap function. 
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Fig. 6. Temperature dependence of the quantity 1/TfT. The inset shows the enlarged part below 
T=10 K. The curves indicate the results of fitting 1/TfT = a x exp{-A C EF/k B T) + (3 x 0.7 
s _1 K _1 . For details, see the text. 
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